THE INCLUSION OF MOISTURE IN A NUMERICAL 
MODEL OF STEADY-STATE FRONTS 



Clifford James Cornelius 



LiDraiy 

Naval Postrraauaie 



Monterey, California 



bcnoot 
i 93940 



NAVAL POSTGRADUATE SCHOOL 

Monterey, California 




THESIS 

THE INCLUSION OF MOISTURE IN A NUMERICAL 
MODEL OF STEADY-STATE FRONTS 

by 

Clifford James Cornelius, Jr. 

Thesis Advisor: R. T. Williams 

March 1974 

Apphovzd fan. puhtic fiaZiaji t; dl&tnAJbation untimite.d. 

T~ / 5*1 / / ^ 



The Inclusion of Moisture in a 



Numerical Model of Steady-State Fronts 



by 



Clifford James Cornelius, Jr. 
Lieutenant, United States Navy 
B.S. in ED., Seton Hall University, 1964 



Submitted in partial fulfillment of the 
requirements for the degree of 



MASTER OF SCIENCE IN METEOROLOGY 



from the 

NAVAL POSTGRADUATE SCHOOL 
March 1974 



she ^ „ 

C 

e. • 



Library 

Naval Postrraauare bcnoo> 
Monterey, Calitorma 93940 



ABSTRACT 

The numerical f rontogenesis model of Williams (1973) 
is modified to include moisture with its subsequent con- 
densation and release of latent hearting. The turbulent 
diffusions of momentum, heat and moisture are repre- 
sented with various coefficients. The numerical solutions 
show realistic quasi-steady fronts forming within one to 
two days. These solutions are examined and compared over 
a range of the various coefficients, and various ranges of 
temperature . 

Inclusion of moisture in the model causes intensifi- 
cation of baroc linicity at mid and upper levels. Also 
noted is a sensitivity of the moist model to reference 
potential temperature due to the exponential relationship 
between saturation vapor pressure and temperature. 



2 



TABLE OF CONTENTS 



I. INTRODUCTION 8 

II. BASIC EQUATIONS AND BOUNDARY CONDITIONS 10 

III. INITIAL CONDITIONS 20 

IV. CONVECTIVE ADJUSTMENT, LATENT HEATING 

AND PRECIPITATION 26 

V. NUMERICAL SOLUTIONS 31 

VI. CONCLUSIONS AND RECOMMENDATIONS 44 

LIST OF REFERENCES 72 

INITIAL DISTRIBUTION LIST 73 

FORM DD 1473 79 



3 



LIST OF TABLES 



I. Listing of Experiments ------------------ 46 

II. Synopsis of Experimental Results ------------- 47 



4 



LIST OF ILLUSTRATIONS 



1. Time variation (days) of the frontal scale (km) at 

p = 99.1 cb for experiments 1 and 2------------ 48 

2. Frontal slope as p(cb) versus y (km) for experiments 

1 and 2 f 49 

3. Time variation (days) of the frontal scale (km) at 

p = 99.1 cb for experiments 12 and 13 - -- -- -- -- --50 

4. Frontal slope as p(cb) versus y (km) for experiments 

12 and 13 51 

5. Cross section of u (m/sec) for experiment 1-.- - -- -- --52 

6. Cross section of u (m/sec) for experiment 2- - -- -- -- - 52 

7. Cross section of 0 , (°K) for experiment 1- - -- -- -- -- 53 

8. Cross section of 0 1 (°K) for experiment 2- - -- -- -- -- 53 

9. Cross section of q(g/kg) for experiment 1---------54 

10. Cross section of q(g/kg) for experiment 2--------- 54 

-3 

11. Cross section of w (xlO mb/sec) for experiment 1- - - - - - 55 

-3 __ 

12. Cross section of w(xl0 mb/sec) for experiment 2- - - - - - 55 

13. Cross section of u (m/sec) for experiment 4- - -- -- -- - 56 

14. Cross section of 0'(°K) for experiment 4- - -- -- -- -- 56 

15. Cross section of q(g/kg) for experiment 4--------- 57 

16. Cross section of w(xl0 mb/sec) for experiment 4- - - - - - 57 

17. Cross section of u (m/sec) for experiment 5- - -- -- -- - 58 

18. Cross section of u(m/sec) for experiment 6- - -- -- -- - 58 

19. Cross section of 0 1 (°K) for experiment 5- - -- -- -- -- 59 

20. Cross section of 0’ (°K) for experiment 6- - -- -- -- -- 59 

21. Cross section of q(g/kg) for experiment 5--------- 60 



5 



22. Cross section of q(g/kg) for experiment 6--------- 60 

23. Cross section of w(xl0 ^mb/sec) for experiment 5- ----- 61 

24. Cross section of w(xl0 ^mb/sec) for experiment 6- ----- 61 

25. Cross section of u (m/sec) for experiment 11 -------- 62 

26. Cross section of 0 , (°K) for experiment 11 --------- 62 

27. Cross section of q(g/kg) for experiment 11- --------63 

-3 

28. Cross section of w(xl0 mb/sec) for experiment 11 ----- 63 

29. Cross section of u (m/sec) for experiment 12 -------- 64 

30. Cross section of u (m/sec) for experiment 13 -------- 64 

31. Cross section of 0' (°K) for experiment 12 --------- 65 

32. Cross section of 0' (°K) for experiment 13 --------- 65 

33. Cross section of q(g/kg) for experiment 12- -------- 66 

34. Cross section of q(g/kg) for experiment 13- --------66 

-3 

35. Cross section of w(xl0 mb/sec) for experiment 12 ----- 67 

36. Cross section of w(xl0 ^mb/sec) for experiment 13 ----- 67 

37. Cross section of u (m/sec) for experiment 14 -------- 68 

38. Cross section of u (m/sec) for experiment 15 -------- 68 

39. Cross section of 0* (°K) for experiment 14 --------- 69 

40. Cross section of 0* (°K) for experiment 15 --------- 69 

41. Cross section of q(g/kg) for experiment 14- --------70 

42. Cross section of q(g/kg) for experiment 15- --------70 

_3 

43. Cross section of w(xl0 mb/sec) for experiment 14 ----- 71 

_3 

44. Cross section of w(xl0 mb/sec) for experiment 15 ----- 71 



6 



ACKNOWLEDGEMENTS 



The author wishes to express his appreciation to 
Dr. Roger Terry Williams for his technical assistance, 
recommendations, encouragement, and for the use of his 
steady-state f rontogenes is model which was modified 
for the experiments conducted. 

Thanks are expressed to Dr. Robert L. Haney, who 
acted as second reader, for his constructive comments on 
the manuscript. 

Thanks are also expressed to Dr. George J. Haltiner 
who provided information on schemes for dry and moist 
convective adjustment and to William Ehrman of the W. R. 
Church Computer Center who unselfishly gave of his time 
to help in improving the author’s programming efficiency. 

The author furthermore wishes to acknowledge the 
W. R. Church Computer Center at the Naval Postgraduate 
School on whose computer system IBM 360/67 all computa- 
tions were performed. 

A final acknowledgement is given to Charlotte S. 
Cornelius who now awaits the author's return after a 
long absence. 



7 



I 



INTRODUCTION 



Hoskins and Bretherton (1972) have shown analytically, 
and Williams (1967, 1972) has shown numerically that dis- 

continuous fronts can form within a finite period of time 
if no turbulent -diffusion is present. These studies sug- 
gest that a discontinuity will form within 24 to 36 hours 
when reasonable initial conditions are used. If turbulent 
diffusion is present, it can be expected that a balance 
will be achieved between the f rontogenetic advections and 
turbulent diffusions of momentum, heat, and moisture. The 
front should remain in this state of quasi-balance as long 
as the large scale deformation field causes f rontogenetic 
advections around the front. 

The numerical f ron togenesis model of Williams (1973) 
shows the effects of horizontal and vertical turbulent 
diffusions of momentum and heat in the formation of these 
quasi-steady fronts. This model, which also includes an 
Ekman boundary layer, will hereafter be referred to as 
W73 and the model as it existed in Williams (1972) will be 
referred to as W72. 



This investigat 


ion 


appends 


moisture with its 


subse- 


quent 


condensation 


and 


release 


of latent heating 


to W73. 


It is 


assumed that 


all 


the late 


nt heat released 


in this 



condensation is used to heat the air isobar ical ly . 
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The purpose of this st 
frontal solutions from the 
period. The physical mode 
one of the models treated 
and Bretherton (1972). Th 
numerically describes fron 
nondivergent horizontal wi 
ing deformation. This def 
in time and is independent 
boundary layer. The moist 
with diffusion are used in 
tities are functions of Y 
from the model is in the f 
in the Y-P plane. 

W73 has been transform 
obtain more accurate value 
equations that have been i 

In Section 2, the basi 
oped, transformed to press 
cation which keeps the pro 
The initial conditions are 
pressure coordinates in Se 
processes and latent heat 
in Section 4. The experim 
Section 5 and conclusions 



udy is to obtain and examine 
model throughout a four day 
1 is essentially the same as 
by Hoskins (1971) and Hoskins 
e basic model, as it now exists, 
togenesis which is forced by a 
nd field which contains stretch- 
ormation wind field is constant 
of height, except in the Ekman 
hydrostatic primitive equations 
which the time dependent quan- 
and P only. Therefore, output 
orm of vertical cross-sections 

ed to pressure coordinates to 
s of pressure for the moisture 
ncluded in the present model, 
c forecast equations are devel- 
ure coordinates and a simplifi- 
blem two dimensional is discussed. 

presented and transformed to 
ction 3. Convective adjustment 
release in the model are included 
ents undertaken are examined in 
are given in Section 6. 
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II 



BASIC EQUATIONS AND BOUNDARY CONDITIONS 



The hydrostatic 
and moisture may be 



primitive equations with diffusion 
written as follows: 



3V 3 (wv) 

+7 ‘ (-V* Sf~ + 

|| +V (0v) + 

|3 ( q v) + |j(*q)- 



3 2 v 

V(j> + f (KxV)=A m V 2 V+C m J , 

ds 



V 20+C e 0 + V 2 « + H c 



A V 2 q+C + M - M / 

q q 3s 2 A ° 



»•! ♦. It ■ 0 < 

|| - L ( s ) 0 , 



( 2 . 1 ) 

( 2 . 2 ) 

(2.3) 

(2.4) 

(2.5) 



where the vertical coordinate has been replaced by s to 
aid in later discussions on transformation. The equations 
are from W73 with slight modification due to the addition 
of the equation of conservation of water wapor (2.3) which 
is analogous to the thermodynamic equation (2.2). The 
function L(s) in the hydrostatic equation (2.5) is a con- 
stant g/0 in height coordinates. Also sinks and sources 
o 

of heating and moisture have been added to the thermody- 
namic equation (2.2) and the conservation of water vapor 
equation (2.3) . 

The notation used is basically the same as W72 which 
uses height as the vertical coordinate. However, in our 
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subsequent transformation to pressure coordinates, 0 is 
defined as a departure from the domain average of the 
initial state, <{> represents geopotential, gz and w is 
ds/dt in the continuity equation (2.4), thereby allowing 
compressibility . 

The quantities, A , A n , and A are the horizontal tur- 

m U q 

bulent diffusion coefficients of momentum, heat, and 

moisture, respectively, while C , C^, and C are the cor- 

m 0 q 

responding vertical coefficients. These constant coeffi- 
cients were assigned various reasonable values in the 
different moisture experiments. Dry and moist convective 
adjustment processes are represented by the functions 
and while latent heating due to large scale condensa- 
tion is represented by H . In the moisture equation, M 

L A 

is a sink or source of moisture due to moist convective 
adjustment while is a sink for moisture due to large 

scale condensation. 

Boundary conditions for the above equations are such 
that the domain is bounded by two rigid quasi-horizontal 
planes. The no-slip condition at the lower boundary and 
the no-stress condition at the top roughly represents the 
troposphere. The conditions are as follows: 



w(x,y,o,t) = w(x,y,H,t) 
V (x ,y,o,t) = 0 



= 0 



3 v 

Js 

36 



c m aI (x 'y' H 't) 



= 0 



c e it (x ' y ' 0 ' t) 

C q |f (X ' y ' 0 ' t) 



30 



C 0 ^ (X ' y ' H 't) 

C q 



= 0 



= 0 



, ( 2 . 6 ) 
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where H represents the distance between the upper and 
lower boundaries. The lower boundary conditions are realis- 

i 

tic for fronts over level terrain. Therefore, fronto- 
genesis in the lower troposphere is more realistic in the' 
model than are conditions near the "tropopause" . The addi- 
tional boundary conditions for moisture are completely 
analogous to the heating conditions in that fluxes of the 
quantities do not pass through the boundaries. 

The approximate steady state solutions to the hydro- 
static primitive equations from [ (2.1) — (2*5) ] are: 

V = U E D{[x(l-e Bs cos Bs)+ ye Bs sin Bsl i 

- [y(l— e Bs cos Bs)- xe Bs sin Bs] j} 

w = 0 , 

2 2 2 

= $ = -D (x + y ) / 2 - f D xy + F(s) 

0 = G (s) 
q = H (s ) 

1/2 

where B = (f/2C^) , D and f are constant and where F(s) 

is such that the hydrostatic equation (2.5) is satisfied. 

The functions G(s) and H(s) are arbitrary functions of s 
since w = 0. These solutions satisfy the boundary condi- 
tions (2.6) and allow for inclusion of the Ekman boundary 
layer. These relations (2.7) satisfy the steady-state equa- 
tions if some small advection terms in the boundary layer 
are neglected. These small terms are usually dropped in 
the Ekman Theory. 
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The equations and the model become two dimensional 
after a subdivision of the dependent variables as follows 



V = U(x,y,x) + u(y,s,t)i + v(y,s 
w = w(y,s,t) 

0 = 0 (y , s , t ) 

<f> = $ (x,y) + TT (y,s,t) 

q = q(y#s, t) 




( 2 . 8 ) 



It is noted that all departures from the approximate 
steady deformation solutions (2.7) are assumed to be in- 
dependent of x. If we substitute the relations (2.8) in 
the i component of (2.1) we obtain: 



3u ^ 3 , . ^ 3 , , ^ rw , -Bs ^ , 

— (uv) + t-(wu) + uD ( 1-e cos Bs ) 

3t 3y ds 

r -Bs _ . -Bs . ^ . 3u ^ -Bs . 

-D[y(l-e cos Bs)-xe sm Bs]t- + vDe sin Bs 

dy 



3 — B s — B s 

+ wD -r — [xe (1-cos Bs)+ye sin Bs]=fv+A 

o s m 



3 2 u 

3y 2 



+ c 



3 2 



u 



m 3s 2 



(2.9) 



Note that the assumption that u is independent of x 
for all time is violated in (2.9) in the two terms where 
x appears. This indicates that the other quantities 
which were assumed independent of x are, in fact, x- 
dependent and this would bring other terms into the equa- 
tionSo In this study we will neglect these added terms 
and apply (2.9) and the other equations at x = 0. If the 
initial u,v,w,0,7T and q fields are independent of x, the 
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error from this approximation will grow slowly and will 
be confined to the boundary layer since the x-dependent 
terms in (2.9) are zero outside the boundary layer. 

Since f r ont ogene s i s occurs very rapidly it is expected 
that the development of x-variations in the dependent 
variables would have only a small effect on the resulting 
quasi-steady front. In any case these effects could not 
be observed in the atmosphere, because atmospheric fronts 
have some variation in the basic fields along the front 
which would be much more important. 

When equation (2.9) is evaluated at x = 0 it becomes 



9u 9 (uv) 9 (wu) _ 9_V f 9u 

3t 9y 3s ~ l u 3y 3y 



v3 V . -Bs . 

— (e sin 

3y 



Bs ) 



8 — B s 

-Vw (e sin Bs)=fv + A 
d s m 



8 2 u 



8y 



+ C 
2 m 



3 2 u 
3s 2 



( 2 . 10 ) 



g s 

where V = -Dy and T = 1-e cos Bs • The j component of 
(2.1) applied at x = 0 takes the form 



3 v 3 . , , 3 , . „3 . 3V -Bs 

_ + w (vv) + _(wv) + rg 7 (W)-Ug 7 e sin 

- (wV)-| — (e B cos B )= - -fu+A — — — • +C 
d s dy m ^ 2 m 



Bs 

8 2 v 
8s 2 



( 2 . 11 ) 



In a similar manner from equation (2.2), the following 
equation for the departure of potential temperature can be 
derived : 

39 . 3 ( v0 ) . 3 (w0 ) . p „ 30 , 9 2 0 , „ 3 2 0 

3t 3y 3s 3y 0 9y 2 0 9g 2 

+ Q + Q + H , (2.12) 

a m c 
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and the following equation for q follows from (2.3), 



9q 3 (vq) 

at a y 




+ 





c 

q 



ill 

as 2 



+ M 



A 




(2.13) 



Equations (2.4) and (2.5) become 



3 v 3 w 

a y + as 



o 



and 



a 7T 

as 



l(s) e 



(2.14) 



(2.15) 



The boundary conditions at the bottom and top of the domain 
respectively are 



= V = 


CD 

u 

II 

> 


_3_0 
3 s 


C H - 
q 3 s 


O 

w 

II 

o 




9u _ 
m 3 s 


c — 

m 3 s 


= w = 


c — = 
e 3s 


3q ^ 

C = 0 , s = h 

q 3s * 


I (2.16) 



If we define the vertical average of a quantity as 

i H 

<()>.= -/ ( ) ds , (2.17) 

ri 

O 

integrate the hydrostatic equation (2.15) with respect to 
s, and remove the vertical mean, it follows that 
. s s 

77 — < 7T > = [ / L ( s ) 0ds- < f L ( s ) 0ds > ] . (2.18) 

o o 

Take the vertical average of (2.14) and use the boundary 
(2.16) which gives 

!_ <v> - 0 . (2.19) 
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This equation states that the total mass flux in the j 
direction is independent of y. In W72 a symmetry argument 
was used to show that this flux must be zero. This argu- 
ment does not hold strictly in the present experiment, 
but it can be expected that at large distances from the 
frontal zone the disturbance mass flux will vanish. Thus, 
it is assumed that 



< v > = o . 



( 2 . 20 ) 



If the development of W72 is used, equation (2.11) can 
be written 



3_v 
3 1 



+ 




— < v v> ] + 




+ 



|-(rvv - <rv> v) 

3y 



3V . -Bs . „ ^ -Bs . „ . „ r 3 , -Bs 

- ^[ue sin Bs-<ue sin Bs> -V [w -^-(e cos 



Bs) 



- <w 



3 , -Bs 

r — (e cos 

3s 



BS) >] 



g 

■r — ( 7 T — < 7 T> ) -f(u-<U>) 

3y 



+ A 



3 2 v 



m 



3y 



+ C [ 
2 m 



3 2 



3 s 



Z + I'F 1 ’ o' 
2 H 3s s=0 



( 2 . 21 ) 



The equations (2.10), (2.12), (2.13), (2.14), (2.18) 

and (2.21) form a complete set which can be solved by a 
pure marching process. The finite difference equations 
conserve mean squares in the advection terms and are de- 
scribed by Williams (1967). The Matsuno scheme is used 
everywhere except in the arrangement of variables and the 
finite difference approximations are the same as those 
used by Williams (1967). In order to close the problem, 
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computational boundaries must be introduced in y. Since 
the disturbance velocities should die out at a sufficient 
distance from the axis of dilation, then 



v( + Y,s,t) = 0 . 



( 2 . 22 ) 



However, there is appreciable inflow across these com- 
putational boundaries since V (+_ Y ) = DY . The quantities 

u, 0 P and q which are advected across the boundaries must 
be specified independent of the interior values if compu- 
tational stability is to be maintained (Platzman, 1954). 
Thus, the following boundary conditions are used: 



u[+(Y + Ay/2),s,t] = u [+(Y + Ay/2),s,0] 

0t+(Y + Ay/2 ) ; s , t] = 0 [+(Y + Ay/2),s,0] 

q[+(Y + Ay/2 ) , s , t ] = q [+(Y + Ay/2),s,0] 



(2.23) 



The computational boundaries y = + Y are placed between 
gridpoints so that the above conditions are actually 
applied at y p +_ (Y + Ay/2). 

The boundary conditions were found to be satisfactory 
except that a solution separation developed near the 
boundaries. This was controlled by using a Matsuno time 
differencing scheme for stability. 

In order to have precise values of pressure at grid- 
points for the moisture experiments it was deemed neces- 
sary to transform W72 to pressure coordinates. 
Minimization of changes to finite differencing schemes 
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of the model have been achieved by defining a parameter s 
which is analogous to a field of z m Specifically, 

s = P Q - P , (2.24) 

a field that will increase in the upward direction and 
from which values of pressure can easily be obtained. 
Further, it is seen that 9s = - 9P so that finite differ- 
encing schemes in W72 can be equally well applied to the 
field s . 

Through the expression in s coordinates the momentum 
equations (2.10) and (2.21), the thermodynamic equation 
(2.12) and the moisture equation (2.13) remain as before 
for the model. However, the continuity equation (2.14) 
in the model actually remains in w and describes a nega- 
tive field of omega in that, 

d s 

w = — = -0) , (2.25) 

Equation (2.5), the hydrostatic equation, differs 
significantly due to the transformation. The derivation 
is briefly as follows. The hydrostatic equation expressed 
in pressure coordinates and the equation of state gives 



i£ = = _ *1 

3p P 



(2.26) 



If a form of the potential temperature definition is 



used, 




(2.27) 
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and combined with (2.26) and (2.27) it follows 



94 > _ R 
3p ~ ” P 



p o 



e . 



(2.28) 



It can be shown that 



R 

P 



C^-) K = c l-(E-) 

p p 9p p 



(2.29) 



The latter form simplifies finite difference energy conser- 
vation. Therefore, substituting equation (2.29) into 
equation (2.28) and transforming to s coordinates we obtain 



= _ c — 
9s p 9s 



C(^) K ] e , 

P o 



(2.30) 



a form of the hydrostatic equation. Note that L(s) is not 
constant as in height coordinates since 



L (s ) 



= *(£_ 
P P. 



c 

p 




(2.31) 



It follows that the complete set of equations solved by 
the marching process remains intact but with L(s) defined 
as in (2.31). 

The boundary conditions (2*16) also remain intact 
except that the lower boundary of vertical motion in the 
model is approximated as 

w = - (jl) = 0 , at s = 0 . (2.32) 
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Ill 



INITIAL CONDITIONS 



The initial conditions in height of s coordinates are 
the same as those employed in W73. The initial potential 
temperature field is given by 



0 (y , z , o) 
where 

a 



86 I H 2 

[t (z- — ) ] -a (— ) arctan (sinh a y) 

dz 2 7T 



fir 



4 



HiL (s) 



96 . 

97 “ 



(3.1) 



The quantity dQ^/dz which is constant is the initial static 
stability, 0 being the initial horizontally averaged 
initial theta and "a" is one-half the total horizontal tem- 
perature variation. 

By substituting equation (3*1) into the th.ermal wind 
equation and integrating we obtain the x-component of the 
rotational wind which is 

u(y,z,o) = -- pp (z - ~) sech a y . (3.2) 

o 

The divergent part of the wind or v field is obtained 
from the quas i-geos trophic equations in z coordinates. 
Details are as in Williams (1972). The remaining fields 
of w and 7T are obtained from the continuity (2.14) and 
hydrostatic (2.15) equations. 
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Once again transformation of the above initial fields 
must be made to pressure coordinates. 

To obtain initial conditions for the theta field we 
must first find 0 ^ , the horizontally averaged initial 

theta based on Qj(p o ) * Derivation is as follows. Given 
knowledge of differentials and the definition of geopoten- 
tial it is seen that 



30 



I 



3p 




dz_ _ 1_ 99 I 84> 

3p g 3 z 3p ' 



(3.3) 



where 30^/3z is the prescribed static 'stability at t = 0 . 

Substitute a definition of potential temperature (2.27) 
to obtain 



!!i . t- I * ( E_, K : e . 

3p 3z g p p I 



(3.4) 



Integrate from p^ to some p , 



f I !!i dp = f . dp 

y 0 9p J g 9z p p o 



(3.5) 



which yields 

0, 



In 



0 I (p o ) 



80 p 

5. *-!■ [(E-) K - ( -^) K ] . 

g 3z K p 0 p o 



(3.6) 



where < = R/C . it follows that 
P _ 

c 80 t 

0 T (p) = 0 T (P ) EXP [-£ 3-^(1- (E-) ) ] 
I I o gdz p 

o 



(3.7) 
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Equation (3.7) is representative of term one of equa- 
tion (3.1). Therefore, the initial potential temperature 

i 

field in pressure coordinates is given by 

C 

e T (y,P,o) = 0 (p ) EXP [-£ V— ( 1 - (r~) )1 

I I O adz P 

O 

-a(2/TT) arctan (sinh ay). (3.8) 



where 
a ~ 

H 



f TT 



v 



30. 



L(s m )( - 9 — } 




(3,9) 



In the model a is given a constant median value by 

solving equation (3,10) at p = p = 60 cb and using the 

m 

static stability at this level. Recall also that 

pressure in equation (3.8) can easily be obtained for 
model purposes as p = p^ - s. 

Further it should be recalled that 0 is a departure 
from the domain average of 0^. This domain average is 
found in the model by finding the vertical average of 
equation (3.7), < 0^>. 

The initial x-component of the velocity or initial 
u equation is found using a form of the hydrostatic equa- 
tion (2.26) and the definition of geostrophic wind. 



1 3<t> / _ _ _ . 

u = - l ' < 3 - 10) 

A form of the thermal wind equation is found to be 



9 u _ 1 d _ _1 3 9 

9p f 9y 9p f p p Q 9y # 



(3.11) 
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60 cb 



If we again assume a median pressure level = 

and that at this level u (p ) =0, we integrate 

m 



p p 

/ 3u _ R 30 / l.p 

3p p f 3y J p p Q 



X K J 

) dp 



(3.12) 



m 



m 



Recalling that u(p^) 



= 0 and that K = 



R/C^ the integra- 



tion yields 



. . p 30 p K 

u(y,p,o) = ft- j~ [ (— ) 

J O 



( ^ L)IC] * 
P o 



(3.13) 



From (Williams, 1972) equation (4.3), it is seen that 



A® 

3y 



2aa 

7T 



sech (ay) . 



(3.14) 



Therefore, the initial u field in pressure coordinates is 
given by 

, x C p 2aa r ,p . K P m x ic, , . 

u(y,p,o) = — — t (— ) - (— ) ] sech (ay) . (3.15) 

p o p o 

The initial divergent wind v is determined as in Wil- 
liams (1972) but using the quas i-geos trophic equation in 
pressure coordinates. The stream function is defined as 



and 



v = 3y/3p 
0) =-3i^/3y 



(3.16) 



and equation (4.10) of Williams (1972) takes the form 



3 2 \p 



30 . 



, T , , , I, 3 *P 2d t , ,30 

+ L ( s ) (- ) = —T L(s ) , 

3p 2 m 8p p t. 3y 2 £ 2 " 3y 



(3.17) 
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where L(s) is assumed 



r / \ f R / K 1 

Ms) = [ — ( ) ] 

P P 

m o 



The equation (3.17) is solved for ^ from the initial 
temperature field (3.8) with the technique of Ogura and 
Charney (1962) and v is then determined from (3.16). 

It is noted that these initial conditions do not 
satisfy the boundary conditions on u, v, and 0 at the 
upper and lower boundaries. A period of adjustment will 
be required to form the surface friction layer and some 
oscillations may be observed later. 

The fields of w and 7T are obtained through the con- 
tinuity and hydrostatic equations (2.14) and (2.15). 

Finally it is necessary to define explicitly methods 
used in determining the initial field of specific humidity 
given (y^PfO) from equation (3.8). 

Saturation vapor pressure e^, is first determined at 
each gridpoint through use of an integrated form of the 
Clausius-Clapeyron equation, an equation expressing satura 
tion vapor pressure as a function of temperature. Accord- 
ing to Hess (1959) 



e 

s 



L : 

.611 EXP [— — ( 
K 

V 



1 

T 




(3.18) 



where e g is expressed in centibars, L c is the latent heat 
of condensation, is the gas constant for water vapor, 
and T q is a constant of integration. 
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Knowing saturation vapor pressure, the saturated speci- 
fic humidity was obtained from the approximated relation 
from Haltiner and Martin (1957) 

. 62 2e 

s , ^ ^ % 



Therefore, by having knowledge of potential temperature 
at t = 0 from equation (3.8), can be found at each point 

using equations (3.18) and (3.19). By defining q to be a 
given fraction of q^ in different experiments, atmospheres 
of different constant initial relative humidities were 
tested . 
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IV. CONVECTIVE ADJUSTMENT , LATENT 
HEATING AND PRECIPITATION 



Inclusion of two additional heating terms in the moist 
thermodynamic equation (2.3) necessitated some further 
amendment to W73. Atmospheric heating by dry convective 
adjustment, , was included in W73 to eliminate any in- 
tense grid scale convection which developed due to the for- 
mation of any dry unstable lapse rates. This redistribution 
of heating was accomplished using averaging techniques to 
ensure that the dry adiabatic lapse rate was not exceeded 

Lapse rates were tested in each vertical column of the 

3 0 

model. Upon encountering an unstable layer (^— > 0), a 

dp 

cumulative average was taken and applied to the gridpoints 

3 0 

within this layer until neutral stability (-^ — = 0) was 

dp 

attained. Simultaneously the model was checked for insta- 
bilities developing below the convectively adjusted layer 
and potential temperatures at these levels were also cumu- 
latively averaged and reassigned the new values. The final 
result of the technique is that the heating is redistri- 
buted to simulate physical processes of the atmosphere and 

3 0 

mathematically that 0 throughout each column of the 

model. A completely analogous scheme was used for saturated 
parcels in the moist experiment, only using equivalent 
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potential temperatures, 0 , in place of potential tem- 

s E 

perature, 0* 

With the addition of moisture, consideration of large 
scale condensation and the subsequent release of latent 
heating became important. 

According to Haltiner ( ) , changes in moisture and 

temperature due to large scale condensation were deter- 
mined as follows. Where relative humidity exceeds 100%, 
the adjustments in temperature St and specific humidity 
6q are determined by the following equations: 

q + 6q = q g (T + 6 t, p) , (4.1) 

and 

C 6T = - L 6q . ( 4 .2) 

P c 



These equations require that the final state of the 
air be exactly saturated and that excess moisture be con- 
densed isobar ically , releasing the latent heat to the air. 

Equations (4.1) and (4.2) can be expressed in itera- 
tive form as 

9q s 

Aq = q g (l-r) + ( j ^~) p At ; q = rq s , (4.3) 

and C AT = - LAq . (4.4) 

• P 



Eliminating At from (4.3) and (4.4) leads to 



Aq 



(1-r) 



L 9q s 

I + — ( — 

1 c ( 9t 
p 



(4.5) 
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At constant pressure the Clausius Clapeyron equation 



can be expressed as 




(4.6) 



Substituting (4.6) into (4.5) we obtain 



q (1-r) 

s 



(4.7) 




Now Aq may be calculated from (4.7) and the correspond- 
ing value of At from (4.4). 

The first estimates of the adjusted temperature and 



There values in turn were used to obtain improved ap- 
proximations. Iterations were containued until the values 
of T 1 and q 1 at each level matched with sufficient accu- 
racy one hundred percent relative humidity. The changes 
were then applied to the potential temperature and speci- 
fic humidity fields of the model accounting for the term 
of equation (2.2) and of equation (2.3) 

The model was then investigated for unstable moist 
convection and its associated heat of condensation. A 
convective adjustment scheme similar to the dry one was 
used on the saturated parcels so that moist static stability 



specific humidity are, therefore 



T 1 = T + At and q 1 = q + Aq 



(4.8) 
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was always maintained. A field of equivalent potential 
temperature was formulated according to Holton (1972) as 



0 



sE 



0 EXP 



L q 

cs 

C T 
P 



(4.9) 



and the moist convective adjustment was designed to ensure 
that the moist adiabatic lapse rate was not exceeded 

(30 sE /8p 1 0) • 

The heat release, Q , as a result of moist convective 

m 

adjustment was included through knowledge of the change of 
equivalent potential temperature during the adjustment. 

By taking the natural logarithm of equation (4.9) it 
is found that 



L q 

In 0 = In 0 + -r 2 -^ (4.10) 

sE C T 

P 

In a model using pressure coordinates, a constant pres- 
sure surface may be assumed, so by differentiating both 
sides of equation (4,10) and recalling the definition of 
potential temperature gives 



d0 



sE 





d q. 



L 

c 



C T 
P 




(4.11) 



Solve the Clausius Clapeyron equation (4.6) for dq , 

s 

substitute in equation (4.11), and solve for dT. The 
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solution is found to be 



d0 



dT = 



sE 

2 

, q L Lq 

e { k (1 + £ 

SE T R C T 2 V 

v p 



)1 , 



(4.12) 
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where d0 is the change in 0 needed to adjust ~ — 
sE sE dp 



sE 



= 0 . 



Since the atmosphere remains saturated throughout this 
process of moist convective adjustment, the Clausius- 
Clapeyron equation was again used to obtain the correspond- 
ing dq or actually dq^ . Once again the changes were then 
applied to the potential temperature and specific humidity 
fields of the irjodel accounting for the term of equation 

(2.2) and M of equation (2.3). 
c 

A cumulative amount of precipitation was then found 
using the relation 



P = 



r p ° 

I <5q dp 



(4.13) 



o 

in agreement with the changes in specific humidity due to 
condensation and moist convective adjustment. Precipita- 
tion was expressed in inches of rainfall in the model. 
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V. NUMERICAL SOLUTIONS 

All numerical results to be shown use the following 
values for the constant: 



At 




540 


sec p = 

*o 


100 


cb p = 60 cb 

m 


H = 80 


cb 


Ap 


= 


1.78 


cb 


Y = 


1800 


km Ay = 60 km 


g - 9. 


-2 

81 m sec 


f 


= 


io” 4 


-1 

sec D = 


io' 5 


sec” 1 0 = 3 00 ° K 
o 






e i 




) = 


270, 


280, 


290, 


3 00 ° K C = 1003 

p 


j oules 


k g~ 1 ° K~ 1 


R 

V 


= 


461 


joules kg 


1 ° K~ 1 


R = 287 joules 


kg~ 1 ° K 


-1 



6 -1 

K = R/C L = 2.5 x 10 joules kg T = 273°K 

pc J o 

a =12 . 5 6 ° K RH = 50,70,85% 




8z 



4 ° K km 



(^— =-) = -.53448, .6, . 8° K cb 

8 P Pm 



A =A a =A =3x 10^, 10^ m^sec ^ 

m 0 q 

2 -1 

= 0, .001409, .0001761 cb sec (which corresponds to 

2 -1 

0, 10, 1.25 m sec in height coordinates) 

2 -1 

= 0, .001409, .0001761 cb sec (which corresponds to 

U 

2 -I 

0, 10, 1.25 m sec in height coordinates) 

2-1 2-1 
= 0, cb sec (which corresponds to 0 m sec in height 

coordinates ) 

where more than one indicated value means that the constant 
was varied in different experiments. Whenever possible, 
the values are the same as those used in W73. 

In this section numerical f ron togenes i s solutions with 
moisture are examined at t = 4 days and compared to the 
corresponding dry atmosphere. 
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As can be seen from Table I, 26 experiments were per- 
formed on five cases, cases selected from Williams (1973) • 
Case 1 contains only horizontal turbulent diffusions of 
momentum, heat, and moisture while cases 2 through 5 con- 
tain both horizontal and vertical turbulent diffusions 
of these parameters. Values chosen for these turbulent 
diffusion coefficients were computed to correspond with 
values used in Williams (1973). 

In Table I, is the reference potential tempera- 

ture ( ° K ) as introduced in equation (3.7). The condensa- 
tion interval indicates the number of time steps between 
successive applications of condensation and moist convec- 
tive adjustment to the model. Turbulent diffusion coeffi- 
cients are then given as previously defined followed by 
the initial relative humidity (percent) and the initial 
static stability (°K/cb). 

Experiments were first performed with a reference 
potential temperature (6^ (P Q ) ) °f 300°K to directly compare 
with the dry solutions of W73. This corresponds to a tem- 
perature range of 62-104°F at the surface along a North- 
South axis. Therefore, experiments were also conducted 
with 0j(P Q ) = 270 , 280 , 290°K to simulate physically realis- 
tic situations. Changes to the dry atmosphere with the 
different references were almost negligible, but due to 
the exponential function in the Clausius-Clapeyron equa- 
tion important differences were observed in the moist 
atmosphere . 
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Experiments were also conducted using different values 
for initial relative humidities and initial static stabili- 
ties o 

A brief tabulation of experimental results is given in 
Table II, comparing widths of frontal zone, frontal slopes, 
vertical motions and precipitation. 

A reasonable* measure of the width of the frontal zone 
is given by 



, - 0 (-Y,s,o) - 0 ( Y , s , o ) 

= |Hi 

1 3y 1 max 



(5.1) 



3 0 

where ds approximated by a one-sided difference. This 

width (km) is indicated at three levels for each experiment 
in the table. 

3 0 

The North-South location of 1 -?r — I at three levels is 

1 3y 1 max 

y (km) and is used as a position indicator of the frontal 
zone where | | 1800 km. 

The pressure level and magnitude of maximum vertical 
motion is then given followed by the North-South position 
of maximum precipitation (km), | \ <_ 1800 km, and the 

corresponding amounts of cumulative precipitation over the 
four day period. 

It is noted from Table II that widths of frontal zones 
between corresponding moist and dry experiments are pri- 
marily effected at mid to upper levels with only small 
deviations being noted at the surface. The expected in- 
crease of thermal gradient at these mid to upper levels 
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due to latent heating in the moist cases is seen to be 
the cause. 

Comparisons of frontal slopes of corresponding experi- 
' * 

ments revealed deviations again at mid levels. The slopes 
appear more vertical throughout these levels in the moist 
experiments. Graphical illustrations will later be pre- 
sented. 

Maxima in vertical motion seem to intensify and be 
found at higher levels along the frontal zone in the cor- 
responding moist cases. This effect can also be directly 
tied to releases of latent heating in this region. 

Cumulative precipitation amounts in Table II are ob- 
served to be primarily functions of initial relative 
humidity and reference potential temperature. 

Primary emphasis on analyses of experimental results 
was on contour plots from the IBM 360/67 computer. These 
results obtained for all 26 experiments are analyses at 
t = 4 days of the fields of U/qjB* and w where 0 * = 0 - 0 ^ 

and 0^ is the horizontally averaged initial potential tem- 
perature. The entire vertical plane is shown from 100 to 
20 cb in these plots. However, the outer portion of the 
domain between |y| = 1200 km and |y| = 1800 km is not shown 

Also note that isolines of specific humidity in the figures 
to follow have different intervals dependent upon reference 
potential temperatures in order to allow for proper defini- 
tion of the field. 
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Due to the large number of contour plots obtained, ex- 
periments 1,2,4,5,6,11,12,13,17,18 were chosen as repre- 
sentative samples for detailed discussion. Experiments 
1,2, 4, 5, 6 and 11 contain only horizontal turbulent diffu- 
sions of momentum, heat, and moisture while experiments 
12,13,17 and 18 contain both horizontal and vertical turbu- 
lent diffusions *of these parameters. 

Figure 1 comparing experiments 1 and 2 shows the 
frontal scale at 99.1 cb reaching steady-state at approxi- 
mately t = 2 days. Little deviation is seen between the 
wet and dry experiment. Figure 2 compares the frontal 
slopes of the same two experiments. As mentioned before 
slopes are seen to be largely the same with the only signi- 
ficant deviation being at midlevels where the moist slopes 
approaches the vertical and then recurves in agreement 
with the dry solution. 

Again this deviation in the wet experiment is probably 
due to the release of latent heating at these levels. 
Similar results were achieved for other cases with more 
extreme deviations being observed at higher reference 
potential temperatures . 

Contour plots for experiments 1 and 2 are shown as 
figures 5 - 12. Figures 5 and 6 show distributions of the 
rotational wind, u for a corresponding dry and moist experi 
ment respectively. It is noted that low level features are 
quite similar in the easterly flow. However, horizontal 
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and vertical wind shear is more intense at low mid to mid 



levels in the warm air in the moist experiment. It is note- 
worthy that this property is dominant throughout the moist 
solution . 

Figures 7 and 8 contrast values of 8' throughout the 
domain* Noted is the more intense gradient of potential 
temperature through the mid levels and the pocket of warm 
air at mid levels most likely due to release of latent 
heating in the moist experiment. 

Fields of specific humidity are contrasted for experi- 
ments 1 and 2 in figures 9 and 10* In the dry experiment 
moisture is merely advected and horizontally diffused 
throughout the domain showing an interesting though physi- 
cally unrealistic situation.* The processes of condensation 
and moist convective adjustment are only applied in the 
moist experiments. By comparing the isolines of specific 
humidity, it is seen that the processes of condensation 
and convective adjustment are indeed reducing the moisture 
content of the air and that the gradient of specific 
humidity is stronger along the frontal zone in the dry 
experiment where no condensation is taking place. 

Vertical motions are investigated for the experiments 
in figures 11 and 12. Upward vertical motions are stronger 
and occur at higher levels closer to the frontal zone in 
the warm air mass for the moist experiment. Downward ver- 
tical motion in the cold air mass was also slightly inten- 
sified with its maximum also at a higher level. 
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Figures 13 - 28 compare the dry and moist fields of 
u, q , 0 1 and w at reference potential temperatures of 
290, 300 and 270°K. Figures 13 - 16 show moist fields 

(B^tp^ = 2 9 0 ° K ) • The horizontal and vertical shear of 
u in figure 13 is seen to be stronger and to be extending 
to higher levels than in the previous moist experiment. 
Stronger thermal gradients are seen to extend higher in 
the atmosphere in figure 14. Latent heat release occurs 
at higher levels and to greater extent as is seen in the 
upper level isentropic pattern to the warm air side of 
the frontal zone. 

As is predicted by the Clausius-Clapeyron equation 
(3.18), the moisture values (figure 15) are significantly 
greater than in the previous experiments. Again note 
the decreased gradient of moisture along the frontal 
zone in the moist experiment. 

Vertical motions, shown in figure 16, are also somewhat 
intensified, occur at higher levels and have slightly 
more intense gradients along the discontinuity in the warm 
air. The downward vertical motions in the cold air are 
also somewhat increased. 

Figures 17 - 24 contrast dry and moist experiments 
with a reference potential temperature of 300°K. Little 
difference is seen between the dry experiment and the 
earlier dry fields based on 280°K except for a northward 
displacement of the fields at the colder reference poten- 
tial temperature. Due to the exponential nature of the 
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moisture parameter, moist fields showed significant changes 
in the following experiments. 

The rotational wind u, as is shown in figures 17 and 
18, has much stronger horizontal and vertical shear in the 
warm air at low mid to upper levels than in previous moist 
experiments and even has an area of easterly winds at upper 
levels in the warm air. The air is much warmer in the 
upper levels due to latent heat release, as indicated by 
figures 19 and 20. Thermal gradients along the frontal 
zone in the warm air are also stronger in the mid and upper 
levels . 

Gradients in the moisture field are greatly intensified 
and specific humidities are in excess of 28 g/kg at low 
levels to the South. Recall that the "dry" q field (figure 
21) is distributed due to advective and horizontal diffusive 
processes only. Perturbations in the moist experiment 
(figure 22) are probably due to the larger quantities of 
latent heat release observed in this experiment. 

- 3 

Vertical motions (figures 23 and 24) exceed 4 x 10 
mb/sec in the warm moist air. These features occur at 
upper levels and a cell containing downward motion is 
observed at these upper levels to the South and seems to be 
somewhat due to the impermeable upper boundary. 

Fields with a reference potential temperature of 270°K 
are finally examined in figures 25 - 28. Due to minimal 

moisture content at these lower temperatures the fields 
are quite comparable to dry experiments previously shown. 
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Shear in the u field (figure 25) and thermal character- 
istics (figure 26) are also quite similar to dry experi- 
ments, Specific humidity (figure 27) reaches a maximum 
of 5 g/kg allowing for only a small possible latent heat 
release. Finally, vertical motions (figure 28) occur 
mainly at the lower levels and are mere closely tied to 
dry experiments. 

Further experiments with case 1 are included in 
Table I. Due to the repetitive nature of the contour plots 
obtained, they are not included as figures, but some fur- 
ther discussion is appropriate at this time. 

Experiments 3 and 8 were run with varying initial rela- 
tive humidities. Higher initial humidities do lead to 
more latent heat release with its increased thermal gradi- 
ents, wind shear, and vertical motion in the warm air but 
did not seem to warrant further experimentation. 

Experiments 9 and 10 were run with varying static sta- 
bilities but due to extensive similarities in the solu- 
tions, further experiments were not required. 

The condensation interval was generally five time steps. 
The iterative procedure of returning relative humidities 
to 100% worked equally well at every ten time steps, never 
exceeding two iterations. However, five time steps were 
finally chosen to compromise computer time and physical 
reality. Only minor changes were observed in the contour 
plots of u, q, 0* and w. 
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Experiments that follow contain both horizontal and 
vertical turbulent diffusions of momentum, heat, and 
moisture. As in the horizontal case previously discussed, 
values of the diffusion coefficients were computed to 
correspond with those used in Williams (1973), 

Figure 3 compares the frontal slopes for the corres- 
ponding dry and moist experiments 12 and 13. Once again, 
slopes are almost identical with the exception of at mid- 
levels, where in the moist experiment maximum potential 
temperature gradient again shifts toward the warm air and 
then recurves to agree with the dry experiment. 

Figure 4 for the same experiments shows frontal scale 
at 99.1 cb versus time with steady-state being obtained 
at approximately c. t = 2 days. Oscillations arise from the 
development of the surface boundary layer, which is not 
present in the initial conditions as was described in 
Williams (1973). Very little difference is observed between 
the moist and dry solutions at low levels which reinforces 
the idea that the model remains essentially unchanged at 
these levels with addition of moisture. 

Contour plots are then shown as figures 29 - 36 for 
experiments 12 and 13. Compared to the corresponding dry 
case, horizontal and vertical wind shear is again increased 
especially at low mid to upper mid levels in the warm air 
(figures 29 and 30). The thermal pattern (figures 31 and 
32) shows increased gradient at the low mid to upper levels 
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in the warm air and the area of maximum latent heat release 



is seen in the isentropic pattern. Vertical influences 
for the moist experiment seem somewhat more pronounced for 
a reference potential temperature of 280°K as compared with 
similar experiments containing only horizontal diffusion. 

Figures 33 and 34 compare the fields of specific 
humidity. The changes continue to be a reduction in 
gradient along the frontal zone in the prime area of latent 
heat release in the moist experiment. Little change is 
observed throughout the rest of the domain. 

Vertical motion fields (figures 35 and 36) for these 
same experiments show low level maxima of 1.5 x 10 3 mb/sec 
in the warm air but an additional more intense maximum of 
vertical motion is found at mid and mid upper levels along 
the frontal zone in the moist experiment. 

Case 3 includes a vertical thermal turbulent diffusion 
coefficient eight times as large as in case 2 previously 
discussed. From Williams (1973) a decrease in-frontal tilt 
is expected along with an increase of frontal scale due to 
the change in diffusion coefficients. This can be observed 
by comparing appropriate figures of the dry atmosphere. 
Other changes to the fields will most likely be due to in- 
clusion of moisture. At any rate our comparisons will 
remain as before, a comparison between the corresponding 
dry and moist experiments. 

Figures 37 and 38, comparing experiments 17 and 18, 
show again the increase of horizontal and vertical wind 
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shear at low mid to upper mid levels to slightly greater 
extent than in the previous moist experiment, experiment 
30. 

The isentropic analyses (figures 39 and 40) show the 
tightening of potential temperature gradient at low mid 
to upper levels once again. However, it is now observed 
that the increased thermal diffusion has smoothed somewhat 
the perturbation in the isentropes in the area of maximum 
latent heat release. This added vertical thermal diffusion 
seems to aid in distributing the i so la t ed . hea t ing induced 
"pocket" that had been noted previously,, 

Figures 41 and 42 compare specific humidity fields of 
the same experiments showing the weakening of gradients at 
low mid to upper levels due to condensation. 

The fields of vertical motion, figures 43 and 44, are 
similar in results to the previously discussed experiment 
containing horizontal and vertical turbulent diffusions 
except for some smoothing of features due to the added 
vertical thermal diffusion of the present experiment. Once 
again the maxima noted in the dry experiment 18 is main- 
tained with an added maximum in vertical motion in the mid 
level warm air along the frontal zone in the moist experi- 
ment . 

After an analysis of the before mentioned experiments, 
it becomes obvious that an observational pattern has 
developed and that further discussion using the line of 
attack would be largely repetitive. 
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Careful analysis of the contour plots of the remain- 



ing of the 26 experiments contained essentially the 
features as those discussed herein. 



same 
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VI 



CONCLUSIONS AND RECOMMENDATIONS 



The numerical model of Williams (1973) is modified 
in this study to include moisture. The numerical solutions 
become quasi-steady within 1-2 days as in W73. Intensifi- 
cation of the baroclinic zone at mid to upper levels seems 
to be the significant addition to the model due to the in- 
clusion of moisture. This intensification of the baro- 
clinic zone at the mid to upper levels is consistently 
manifested by the increases of thermal gradient, wind 
shear, and vertical motions in the warm air along the 
frontal zone which are due to latent heat release. Second- 
ly, it should be noted that the model is essentially un- 
changed at low levels with the inclusion of moisture. 
Thirdly, the model with moisture shows dependence upon 
the reference potential temperature, chosen due 

to the exponential relationship between moisture and tem- 
perature as is indicated by the Claus ius - Clapeyr on equa- 
tion (3.18). It is recalled that the dry model is 
relatively insensitive to choice of this reference tem- 
perature. The moist model illustrates this exponential 
thermal dependence in the various experiments at the 
reference potential temperatures of 270, 280, 290, and 

3 0 0 ° K . 

In future studies using Williams model, attention 
should be directed towards transformation to In p 
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coordinates for a better depiction of upper level gradi- 
ents, Upper level boundaries might also be changed and 
raised to include a stratosphere and finally a more 
sophisticated convective scheme might be applied to the 
model such as that of Arakawa and Schubert (1973). 
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TABLE I. Listing of Experiments 
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Figure 1. Frontal scale. Experiments 1 and 2. 
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Figure 2. Frontal slope. Experiments 1 and 2, 
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Figure 3. Frontal scale. Experiments 12 and 13. 
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Figure 4. Frontal slope. Experiments 12 and 13. 
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Figure 5. u(m/sec) = Exp. 1. 
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Figure 8. 0 ' (°K) . 
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Figure 16. w(xlO ^mb/sec) . Exp. 4. 
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Figure 18. u(m/sec). Exp. 6, 
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Figure 21. q(gAg). Exp. 5. Y(l<m) 




Figure 22 0 q(g/kg). Exp, 6* Y( Icm) 
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Figure 25. u(m/sec). Exp. 11. 
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Figure 27. q(g/kg). Exp. 11. Y(km) 
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Figure 28. w(xlO mb/sec). Exp. 11. 
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Figure 29. u(m/sec). Exp. 12. Y(km) 




Figure 30. u (m/sec). Exp. 13. Y(km) 
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Figure 33. q(gAg)* Exp. 12„ Y(km) 
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Figure 34 . 
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Figure 35. w(xlO mb/sec). Exp. 12, 
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Figure 36. w(xl0 mb/sec) « Exp. 13. 



67 



P(cb) p(cW 




Figure 37 . u(m/sec). Exp. 17. 



Y{km) 



20 • 



40 - 



60 * 



80 * 



100 * 



Cl trx-|0 





12 00 900 600 300 

Figure 38 0 u(m/sec) . Exp. 18« 



0 -300 

Y(km) 



i r i 

-600 -900 -1200 



i 



68 



p(cb) * P(C« 




Figure 39. 0'(°K). Exp. 17. 
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Figure 42. 
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Figure 43. w(xlO mb/sec). Exp. 17. 
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